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ABSTRACT 

We have recently identified the widest very low-mass binary (2M0126AB), 
consisting of an M6.5V and an M8V dwarf with a separation of ~5100 AU, 
which is twice as large as that of the second widest known system and an order 
of magnitude larger than those of all other previously known wide very low-mass 
binaries. If this binary belongs to the field population, its constituents would 
have masses of ~0.09 Mq, at the lower end of the stellar regime. However, in the 
discovery paper we pointed out that its proper motion and position in the sky are 
both consistent with being a member of the young (30 Myr) Tucana/Horologium 
association, raising the possibility that the binary is a pair of ~0.02 Mq brown 
dwarfs. We obtained optical spectroscopy at the Gemini South Observatory 
in order to constrain the age of the pair and clarify its nature. The absence 
of lithium absorption at 671 nm, modest Ha emission, and the strength of the 
gravity-sensitive Na doublet at 818 nm all point toward an age of at least 200 Myr, 
ruling out the possibility that the binary is a member of Tucana/Horologium. 
We further estimate that the binary is younger than 2 Gyr based on its expected 
lifetime in the galactic disk. 

Subject headings: Stars: low-mass, brown dwarfs — stars: individual (2MASS 
J012655.49-502238.8, 2MASS012702.83-502321.1) 
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Introduction 



Binarity is ubiquitous in stellar systems, from the most massive stars to the substellar 
regime. The statistical properties of binaries and higher multiple systems retain information 
on the physical processes that led to their formation long after they have dispersed beyond 
their star-forming regions. While n early all early-B sta.rs res ide in binary systems, with a 
binarity fr action as high as ~80 % (IKouwenhoven et al.l 120051 ). this fraction falls to as low 
as ~20% (IBurgasser et al.l 120071 ) for brown dwarfs. Also, while stellar binaries are found 



with separ ations reaching 35 000 AU for early-B primaries and 20 000 AU for A primaries 
( lAbtlll988l ). only a handful of substellar systems have separations beyond 200 AU. Most of 
these have been found in star-forming regions and open clusters, such as Op h 1623-2402 



(212 AU) and Oph 1622-2405 (243 AU) in the star-forming clouds of Qphiuc hus fl Close et a 



20071 ). 2M1101-7732 (241.9 AU) in the Chamaeleon I star-forming region flLuhmanl l2004l ) 



and the system SE70/ S Ori68, a likely 1700 ± 300 AU-wide system in the a Orionis cluster 
( jCaballero et al.ll2006l). Only three wide very low mass (VLM) binaries are known in the 
field, DENIS055- 44 jBil leres et al.l I2OO5I ) (~220 A U), K5nigstuhl 1 A B (1800 ± 170 AU) 



(1Caballeroll2007bl ) and 2M0126AB (5100 ±400 AU) flArtigau et al.ll2007l ). ICaballerol fl2007al ) 
established that wide late-type binaries with a mass ratio >0.5 are rare objects in the field, 
with only 1.2 ± 0.9% of VLM stars and BDs residing in such systems. 

The 2M0126AB system represents the most e xtreme case of VLM binary known. It 
wa s discovered as a c ommon proper motion pair by lArtigau et al.l (120071 ). and later verified 
by ICaballerol (l2007a[ ). The derived probability of having, over the whole sky, a single pair 
of unrelat ed late-Ms with ma. tching proper motion and distances within our uncertainties 
is ~0.002 (lArtigau et al.ll2007l ): hence the two components are almost certainly gravitation- 
ally bound. Its c omponents have spectral types of M6.5 and M8 as confirmed by GNIRS 
( jElias et al.l l2006l ) spectroscopy. While near-infrared spectroscopy shows that KI equiva- 
lent widths are compatible with those of field objects, th e pair falls near the core of the 
Tucana/Horologium (Tuc/Hor) (IZuckerman fc Songj|2004l ) association and shares its bulk 
motion, suggesting that 2M0126AB may be a member of this 30 Myr old moving group. If 
indeed it is a member of this group, the low temperatures of its components combined with a 
very young age would imply masses at the lower limit of the brown dwarf rea lm. Otherwise, 



i f this pair is a field object (i.e. an age greater than 1 Gyr) the models of IChabrier et al. 
( 120001 ) indicate masses of 0.095 M© and 0.092 M©, close to the lower limit of the stellar 
regime. 

Here we present new optical spectroscopy of both components of 2M0126AB to constrain 
the properties of this odd pair, either as a member of a young association or as a much older 
field object. The wavelength interval selected contains important age or gravity indicators: 
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the lithium feature at 670.9 nm, Ha and the 820 nm Na doublet. Observations and data 
reduction are described in §[2j Results and a discussion on the physical properties constrained 
by these observations are given in § [Hi 

2. Observations and data reduction 

The dataset described here was obtained with the GMOS-S spectrograph at the Gemini 
South telescope on 2007 October 3. A 0.5"-wide slit was used with the R400 grating and 
OG515 filter; the resulting resolving power was R~1900 for a 600 — 950 nm wavelength 
coverage. The slit was aligned as to obtain a spectrum of both objects simultaneously and the 
observations were made using the nod-and-shuffle mode for improved sky-line subtraction. 
Two 44.7-minute and one 12.2-minute exposures were obtained for a total integration time 
of 101 minutes. The observations were obtained at a mean airmass of 1.5 under photometric 
conditions. The exposures were taken at slightly different grating angles providing a spectral 
dithering in order to fill the 3 nm wavelength-coverage gaps caused by the spacing between 
the three GMOS detectors. The 790 nm, 800 nm and 810 nm central wavelengths for each 
exposure displaced these gap in the 716 — 740 nm and 870 — 882 nm intervals, i.e. in regions 
devoided of diagnostic spectral features. A white dwarf (EG131) was observed on 2007 
September 11 with the same setup for instrumental and telluric corrections. 

The individual nod-and-shuffle frames were first dark and bias subtracted. Each nod- 
and-shuffle frame was then calibrated using flats taken immediately before (790 nm setup) 
or after (800 nm and 810 nm setup) the science observation. The two nod-and-shuffle frames 
were then pair-subtracted, and corrected for trace distortion. Wavelength calibration was 
performed using copper-argon lamp spectra taken as part of the night-time calibrations. 
Spectra were extracted from both the positive and negative nod-and-shuffle traces for each 
component, resulting in a total of 6 spectra that were interpolated on a finer common- 
wavelength grid and median combined. The final signal-to-noise ratio per resolution element 
ranges from about 15 around 650 nm to 70 around 800 nm for 2M0126A. The signal-to-noise 
obtained for 2M0126B is ~20% lower. 

An image of the 2M0126AB field was taken prior to the spectroscopic observation. 
This 30.5-s i-band image used the GMOS-S central CCD without pixel binning, providing a 
2.5' X 5.6' FOV with a 0.073" pixel sampling. The image was bias subtracted, fringe corrected 
and fiat-fielded. This image was used to search for additional comoving objects and test the 
possibility that one or both components of the system is a tight binary itself. 
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Results & Discussion 



Figure [T] shows the spectra of 2M0126A and 2M0126B. Both spectra show the hallmark 
absorption features of late-M dwarfs such as TiO , KI a nd Nal. A comparison with template 
M6.5 and M7.5 spectra from iKirkpatrick et al.l (119991 ) shows little difference with typical 
field objects, the differences around 760 nm and between 930 nm and 950 nm being due to 
telluric absorption (corrected in our spectra, unaccounted for in the reference spectra). 



We derive an optical spectral type using the iMartfn et al.l (119991 ) polynomial relations 
for the PCS index for both component s. These spectral types agree within 0.2 subclass to 
those measured in lArtigau et al.l ( 120071 ) and therefore we keep spectral types of M6.5 ± 0.5 
and M8.0 ± 0.5 for our analysis. Both 2M0126A and 2M0126B show moderate Ha emission 
(see inset in FigureH]) with pseudo-equivalent widths (pEW) of —3.44 A and —7.32 A, respec- 
tiv ely. This measur e ment was converted into log (LHa/i^boi) using the procedure described 
by lWalkowicz et al.l (120041 ): we obtained log (LH«/i^boi) = —4.5 and —4.4 for components A 
and B, respectively. The measured Ha pEW for both components are typical for objects of 
these spectral types (IWest et al.ll2004l ). From these values alone we cannot rule out the pos- 
sibility that the pair is a member of the Tuc/Hor association as young late- Ms have widely 
varying levels of activity; the level observed for 2M0126A is comparable to the least active 
members of the 50 Myr old open cluster IC 2391 while 2M0126B h a s Ha strength similar 
to that of other late-Ms in the sample of iBarrado y Navascues et al.l (120041 ). Measurements 
of the Ha pEW therefore provides little help in establishing the true nature of 2M0126AB, 
excluding neither of the hypotheses. 

The absence of lithium absorption gives the best clue that 2M0126AB is not a member of 
the Tuc/Hor association (see inset in Fig ure [J). At a bolometric luminosity of log (L/Lp,) ~ 
-3.1, both objects should take 200 Myr Jp'Antona fc Mazzitellilll997l : Ichabrier et allboooh 
to destroy lithium, impl ying a minimum ag e that i s nea rly an order of magnitude greater than 
the 30 Myr of Tuc/Hor. IKirkpatrick et al.l (120061 . |2008| ) mention that in some young objects, 
the lithium signature can be absent due to very low gravity and the intrinsically weaker alkali 
lines seen in these objects. This would explain the seemingly contradictory spectroscopic 
features of 2MASS J01415823-4633574: no detectable lithium at 671 nm despite numerous 
low gravit y indicators pointing t oward an age of 5-10 Myr and an estimate gravity of log (g) ~ 
4.0 ± 0.5. IChabrier et al.l (|2000l ) models predict, for objects at the luminosities of 2M0126A 
and 2M0126B, log (g) = 4.2 at an age of 30 Myr (Tuc/Hor scenario), log^f = 5.0 at 200 Myr 
(minimum age with lithium depletion) and log (g) = 5.3 at 1 Gyr (half-life of the binary due 
to tidal disruption in the galactic disk). While the log (g) predicted for 2M0126A and B in the 
Tuc/Hor scenario is close to that of 2MASS J01415823-4633574, which may cast a doubt on 
the validity of the lithium test in this case, we note that objects later than M5 in the 50 Myr- 
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old o pen cluster IC2391 do show significant lithium absorption (IBarrado y Navascues et al. 
20041 ) and that, at the same temperature as 2M0126A and B, IC 2391 late-Ms have surface 
gravities within about 0.4 dex of that predicted for 2M0126A and 2M0126B at 30 Myr. 



The Na I doublet at 818 nm is a good gravity indicator (IKirkpatrick et al.lll99ll . l2006l ). 
The Nal doublet equivalent widths of 7.04 ± 0.12 A and 6.69 ± 0.20 A for 2M0126A and B 
respectively, are consistent with th ose of field objects but about 2 A larger than the values 
found for late-M Pleiades members (IMartm et al.lll996[). This indicates tha t both objects are 
significantly older than the Pleiades (130 Myr; iBarrado y Navascues et al.l (120041 )). which in 
turn implies surface gravities of log (g) > 5.0. Such a high surface gravity further implies 
that the binary system would clearly show Li absorption should it be present. This confirms 
the lower limit on the age set by the lithium test at 200 Myr. We therefore establish that 
2M0126AB cannot be a member of the Tuc/Hor association and is more likely a field pair. 

Interestingly, this pair has an expected survival time in the galactic disk significantly 
shorter than t he age of the disk itse lf. This timescale has been estimated by scaling to 
2M0126AB the IWeinberg et al.l (119871 ) results for binaries in the solar neighborhood. With 
an half-life of ~1 Gyr, we can estimate at the 68% confidence level (equivalent to la) that 
the pair has a i i age between 0.2 and 1.8 Gyr. Masses have been determined from the 
Chabrier et al.l (l2000l ) models for ages of 0.2 Gyr and >1 Gyr and are given in Table [H 



We measured a list of 12 spectroscopic indices (see table [T]) compiled by ICruz fc Reid 



( 12002! ) in order to verify whether 2M0126AB shows any peculiarity compared to field M 
dwarfs. All indicators fall within the trends observed for field objects; three are shown 
in Figure [2l None of the indicators deviate from the overall trends for field M dwarfs. 
In particular, the Ti05, being much stronger in sub-dwarfs and extreme sub-dwarfs, is a 
useful probe of metallicity. Our CaHl, CaH2 and CaH3 versus Ti05 measurements clearly 



put 2M0126A in the non-subdwarf part of Figure 1 in iGizid (119971 ). which extends only 
up to a spectral type of M7. This is in agreement with the <2 Gyr age expected from 
dynamical considerations. Using the relations between t he Ti05 and CaH2 indices versus the 
J-band absolute magnitude given by lCruz &: ReidI (120021 ). we derive a photometric distance of 
82 ± 5 pc for 2M0126A. The photometric distance for 2M0126B cannot be easily determined 
through this method as the trends break for object s later than M7. Th is revised distance 
is somewhat larger than the distance estimated in lArtigau et al.l (120071 ) and would revise 
upward the physical separation of the pair to about 6700 AU. Using the iReid et al.l (119951 ) 
relation the Ti05 index alone points toward a slightly earlier spectral type (M6.0 ± 0.5) for 
2M0126A. This mild inconsistency in the distance to 2M0126AB illustrates the need for a 
more reliable distance measurement based on parallax. 



Using the GMOS i-band image of 2M0126AB, we rule out, using the subtraction of 
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a field stars' PSF, that either of its components is itself an equal-luminosity binary with 
a separation >0.2". An object ~4" North-West of 2M0126A was examined as a plausible 
close-in companion, 2.0 mag fainter than 2M0126A. This object is undetected in J-band 
images taken in 2006 June, implying a flux ratio greater than ~3.0 mag. This object is 
thus at least 1.0 mag bluer in ^ab — J than 2M0126A and therefore cannot realistically be 
a later-type companion. No other object, down to a A^ab ~4.5 mag, is seen within 5" 
of either component. The flux ratio between 2M0126A and 2M0126B is measured to be 
A i = 0.518 ± 0.015 mag, consi stent with, but much more accurate, than the SuperCosmos 



Sky Survey (IHambly et al.ll200ll ) /-band photometry. A separation of 81.89±0.2" is measured 



between the two components, consistent with previous measurements. 

Overall, aside from the fact that they compose the widest VLM binary, 2M0126A and 
2M0126B are basically indistinguishable from the bulk of field M dwarfs, in age, surface 
gravity, activity level or metallicity. The fact that 2M0126AB is a relatively old system 
makes it all the more remarkable as compared to the bulk of wide VLM binaries found in 
star-forming regions. This object demonstrates that at least a handful of dynamically very 
fragile systems do survive to become field objects; a result that future work on the dynamics 
of star-forming regions will have to take into account. 

The properties of 2M0126AB are interresting to compare to those of their more massive 
counterparts. The correlation bet ween the up per limit on separation versus the mass of the 
primary (2500 Mp;.^^^^^) noted by Abt (jl988l ) in the B5-K0 interval clearly doesn't hold for 
the system described here, where we would expect a maximum separation of about 60 AU. 
Interrestingly, the order of magnitude of this limit is consistent with the sharp cutoff at 
~ 50 AU in the discribution of VLM binaries. While the discovery of 2M0126AB shows that 
much wider VLM binaries exist, the next step in the study of such systems would be a census 
of similar systems through a dedicated survey, in which high care is taken in assessing the 
sensitivity and completeness limits such that clear and robust conclusions can be drawn. As 
these systems are likely to be very rare, a large sample of mid-to-late- Ms would be needed. 
It is also noteworthy that among the 500 or so L and T dwarfs in the 2MASS catalog, no 
binary wider than 2" has been identified. How rare are pairs with total masses much below 
that of 2M0126AB at similar separations, if they exist at all, remains to be seen. 

An interesting avenue of investigation will also be to determine whether forming pairs 
like 2M0126AB is fundamentally different from forming VLM binaries separated by tens of 
AUs. One striking property of VLM binaries , as compared to stella r binaries, is the strongly 



peaked mass ratio distribution toward unity (IBurgasser et al.ll2007l ). with 77% of the systems 



having M2/M1 > 0.8, while late-F to K dwarfs binaries show a much flatter distribution. 
Interestingly, 2M0126AB has M2/M1 ~ 0.95, more in line with the VLM binary distribution. 
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While obviously no broad-ranging conclusions can be drawn from the discovery of a single 
object, the discovery of only a handful of systems similar to 2M0126AB should clarify if they 
share the formation mechanisms of their siblings and represent the tail of the distribution, 
possibly dynamically excited during their formation, or if the formation scenarios themselves 
differ. 

Based on observations obtained at the Gemini Observatory (program ID : GS-2007B- 
Q-18), which is operated by the Association of Universities for Research in Astronomy, Inc., 
under a cooperative agreement with the NSF on behalf of the Gemini partnership: the Na- 
tional Science Foundation (United States), the Science and Technology Facilities Council 
(United Kingdom), the National Research Council (Canada), CONICYT (Chile), the Aus- 
tralian Research Council (Austraha), Ministerio da Ciencia e Tecnologia (Brazil) and SECYT 
(Argentina) . DL is supported via a postdoctoral fellowship from the Fonds quebecois de la 
recherche sur la nature et les technologies. RD is financially supported via a grant from the 
National Research and Engenerring Council of Canada. This publication has made use of 
the VLM Binaries Archive maintained by Nick Siegler at http://www.vlmbinaries.org. 
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Fig. 1.— Spectra of 2M0126A (bottom) and 2M0126B (top, offset by +1) as well as two 
template objects, a M6.5 (2MASS J02422+1343, dotted line, overplotted on 2M0126A) and 
an M7.5 (2MASS J2585+1520, dotted line, offset by +1, overplotted on 2M0126B). All 
spectra have been normalized to their median flux over the 810-840 nm spectral interval. 
The inset shows the 650-675 nm interval that contains two important age diagnostic features, 
the undetected lithium line at 670.8 nm and Ha at 656.3 nm. The inset normalization is the 
same as that of the main figure but the spectra have been offset by 0.2 for clarity. 
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Fig. 2. — VO-a, Ti05 and Ha indices as a function of spectral type for the ICruz fc Reid 
(120021 ) sample compared to the values measured for 2M0126A and 2M0126B (data points 
with error bars). The values measured for our binary are consistent with the overall trends 
of the field M dwarf sample, indicating neither particularly large chromospheric activity nor 
signs of peculiar metallicity. 
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Table 1. Parameters of 2MASS J0126AB 
Parameter 2M0126A 2M0126B 



2MASS designation 


J012655.49-502238 


.8 J012702.83-5( 


Angular separation ^ ["] 


81.86 ±0.10 


Ha cos [mas/yr] 


131±9 


135±9 


l^s^ [mas/yr] 


-53±15 


-47±15 


y a 


21.8 


22.2 




17.2±0.3 


17.6±0.3 


AiAB 


0.518 ±0.015 




14.61±0.04 


14.81±0.05 


H" 


14.05±0.05 


14.16±0.04 




13.68±0.05 


13.62±0.05 


FeH 1.20 //m^ [A] 


6.0±1.8 


14.5±1.9 


KI 1.25/xm^ [A] 


9.8±1.1 


10.3±1.0 


Ha pEW [A] 


-3.44±0.40 


-7.32±0.50 


log(LHa/-C/bol) 


-4.5 


-4.4 


Li 671 nm [A] 


<0.27 


<0.52 


Na 1 EW [A] 


7.04±0.12 


6.69±0.20 


CaOH'i 


0.32±0.05 


0.46±0.11 


Ha-i 


1.48±0.03 


2.16±0.13 


CaHl 


0.87±0.07 


0.94±0.11 


CaH2 


0.277±0.007 


0.270±0.012 


CaHS'^'S 


0.581±0.009 


0.563±0.017 


TiO-a<^'^ 


2.71±0.09 


2.81±0.16 


Ti02<i 


0.32±0.02 


0.33±0.03 


Ti03'^ 


0.55±0.03 


0.54±0.04 


7104-^ 


0.53±0.03 


0.71±0.06 


Ti05 -^'S 


0.219±0.008 


0.284±0.016 


VO-a'^'^ 


1.073±0.006 


1.104±0.011 


VO-b^ 


1.256±0.009 


1.361±0.013 


PC3d 


1.569±0.008 


1.850±0.014 


Near-IR^ & optical SpT 


M6.5V±0.5 


M8V±0.5 


Photometric distance ^ [pc] 


63±5 pc 


61±6 
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Table 1 — Continued 



Parameter 


2M0126A 2M0126B 


Physical separation^ [AU] 


5100±400 


Teff - [K] 


2670±180 2490±180 


Mass (200 Myr) ^ [M©] 


0.068±0.005 0.065±0.005 


Mass (> 1 Gyr) [Mq] 


0.095±0.005 0.092±0.005 


Age [Myr] 


>200 


Mbof 


12.76±0.14 12.75±0.14 


Luminosity^ [log(L/L0)] 


-3.21±0.05 -3.20±0.05 



'^From lArtigau et al.l (120071 ). see this reference for more 
details. 

'^From the SuperCosmos Sky Survey catalogue. 
'^From the 2MASS point source catalogue. 



iCruz fc ReidI fl2002h index 



Kirkpatrick et al.l (Il999l ) index 



Martin et all (119991 ) index 



^Gizisl fll997h index 



''Determined using the IChabrier et al.l (120001 ) evolution 
models and the estimated J, H and Kg absolute magni- 
tudes. 



